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摘 要 本 文 对 液 气温 区 Vuilleumier 耦合 脉冲 管 (VM-PT) 制冷 机 进行 了 数值 计算 研究 . 为 了 减 小 来 自 室温 的 漏 热 ， 将 
气 库 置 于 一 级 冷 头 之 上 . 采用 SAGE 软件 建立 了 数值 计算 模型 ， 以 无 负荷 制冷 温度 为 目标 函数 对 制冷 机 运行 参数 进行 了 
优化 . 为 了 进一步 了 解 制冷 机 的 工作 原理 , 研究 了 脉冲 管制 冷 机 内 灼 流 、 声 功 、 炉 流 、 压 力 波 、 体 积 流 率 及 温度 的 沿 程 分 
布 , 并 对 其 相位 关系 和 可 用 能 损失 进行 了 分 析 . 计算 结果 表明 , 该 制冷 机 可 在 平均 压力 0.9 MPa、 压 比 1.61、 运行 频率 0.8 
Hz 时 得 到 3.69 K 的 无 负荷 制冷 温度 , 在 4.2 K 时 可 提供 16.24 mW 制冷 量 . 
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一 Numerical Study of VM-PT Cryocooler Operat t Liquid Helium Temperature 
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Abstract A detailed numerical stu 5 illeumier coupling pulse tube (VM-PT) cryocooler 
operating at liquid helium a t 0 in this paper. In order to reduce the heat leak 
from the room temperature Abts the reservoir was located at the cold end of the first stage. 
A numerical model of the Re was conducted by using the software SAGE. The operating 
parameters were ee. taking the no-load temperature as the target function. In order to 
further study the coupl echanism of the cryocooler, distributions of 


entropy flow, pressure wave, and volume flow rate are presented and 


thalpy flow, acoustic power, 
cussed. Meanwhile, analyzes 
of phase relationship and exergy loss are also performed. A no-Bhd "temperature of 3.69 K could be 
attained with a mean pressure of 0.9 MPa, pressure ratio of 6， operating frequency of 0.8 Hz in the 
simulation and the cooling power at 4.2 K is about 16.24Xkm 。 
Key words liquid helium temperature; VM ms pulse tube cryocooler; numerical simulation 


0 引言 A 
在 低温 超 导 、 低温 电子 学 、 低温 医学 和 空间 技术 


频 运行 条 件 下 , 蓄冷 器 的 流 阻 损失 较 大 , 效率 较 低 。 


等 方面 , 液 氨 温 区 制冷 机 起 着 非常 重要 的 作用 . 目前 
液 氮 温 区 制冷 机 主要 有 G-M 制冷 机 、G-M 型 脉冲 管 
制冷 机 以 及 高 频 斯 特 林 型 脉冲 管制 冷 机 . G-M 制冷 
机 和 G-M 型 脉冲 管制 冷 机 采用 所 气压 缩 机 结合 高 低 
阀 切换 提供 压力 波动 , 不可逆 损 失 导致 其 热效率 较 
低 , 且 需 要 定期 维护 趾 . 同时 , 其 体积 过 大 的 特点 限 
制 了 它 在 空间 领域 的 应 用 . 液 氨 温 区 高 频 斯 特 林 型 
脉冲 管制 冷 机 是 近年 来 的 一 个 研究 热点 ， 目 前 已 采 
用 二 级 到 四 级 结构 获得 了 液 氮 温度 R44， 运行 频率 
一 般 在 20~50 Hz 左右. 液 氨 温 区 制冷 机 一 般 要 采用 


VM 制冷 机 是 斯 特 林 型 制冷 机 的 一 种 ， 和 斯 特 
林 制 冷 机 不 同 的 是 ， 它 利用 热 压 缩 机 进行 驱动 。 理 
论 上 讲 ，VM 制冷 机 具有 和 斯 特 林 制冷 机 相同 的 效 
率 . VM 制冷 机 的 运行 频率 一 般 低 于 5 Hz, 其 可 采用 
较 大 粒 径 的 蓄冷 材料 填充 低温 蓄冷 器 以 减 小 流 阻 损 
失 , 可 实现 较 高 的 蓄冷 器 效率 . 因此 , 研究 液 氮 温 区 
VM 制冷 机 对 发 展 高 效 液 氨 温 区 制冷 机 具有 十 分 重 
要 的 意义 。 

在 10 K 温 区 单 级 VM 制冷 机 的 基础 上 回 ， 潘 
长 钙 等 日 提出 了 VM-PT 型 制冷 机 ， 即 在 单 级 VM 


较 小 的 颗粒 状 或 球状 磁性 蓄冷 材料 填充 蕾 冷 锅 ， 高 
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冷 头 上 气 耦 合 一 级 同 轴 型 脉冲 管制 冷 机 . 为 了 方便 
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调节 , 采用 了 室温 小 孔 - 气 库 进 行 调 相 , 并 最 终 得 到 
了 4.4 KK 的 无 负荷 制冷 温度 . 但 这 种 调 相 方式 会 引 
入 来 自 室温 的 漏 热 ， 恶 化 一 级 冷 头 的 性 能 ， 使 得 一 
级 冷 头 温度 偏 高 ， 进 而 影响 无 负荷 制冷 温度 . 同时 ， 
这 种 制冷 机 的 耦合 机 理 、 相 位 关系 等 也 有 符 于 进 一 
步 研究 . 为 了 获得 更 低 的 制冷 温度 ,本文 开展 了 低 
温 调 相 VM-PT 制冷 机 的 数值 计算 研究 ， 将 气 库 置 
于 一 级 冷 头 之 上 以 减 小 来 自 室温 的 漏 热 . 为 了 对 制 
冷 机 耦合 机 理 和 工作 原理 进行 深入 的 了 解 并 指导 实 
验 系统 的 设计 ， 在 对 制冷 机 运行 参数 进行 优化 的 基 
础 上 , 本 文 开 展 了 能 流 、 相 位 及 可 用 能 损失 分 析 , 并 
对 制冷 机 的 性 能 进行 了 初步 预测 。 


1 制冷 机 绪 松 


图 1 为 低温 调 相 VM-PT 制冷 机 结构 示意 图 . 在 
单 级 VM 制冷 机 的 冷 腔 开 一 小 孔 ， 作 为 同 轴 型 脉冲 
管制 冷 机 的 入 口 , 这 样 , 除了 流动 引起 的 损失 外 , 脉 
冲 管制 冷 机 内 的 压力 和 VM 制冷 机 内 基本 相同 . 为 
了 减 小 来 自 室温 的 漏 热 , 脉冲 管 采 用 低温 -小 孔 气 
的 调 相 方式 . 气 库 采 用 环形 结构 ， 置 于 一 级 ; 
同时 可 作为 脉冲 管制 冷 机 的 冷 屏 . 低温 置 


阀门 ， 因此 小 孔 闪 采用 毛细 管 来 代 殉 oN 
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图 1 低温 调 相 VM-PT 制冷 机 结构 图 
Fig. 1 Schematic of VM-PT cryocooler with cold phase shifter 


表 1 为 制冷 机 初步 设计 参数 ， 后续 数值 计算 也 
按照 此 参数 进行 建 模 . 其 中 , 冷 端 蓄冷 器 采用 两 层 方 
式 进行 填充 , 第 一 层 填 充 45 mm 200# 不 锈 钢丝 网 ， 


第 二 层 填充 80 mm 直径 0.4~0.45 mm 的 铅球 . 脉冲 


管 蓄冷 需 则 全 部 填充 直径 0.2~0.25 mm 的 EraNi, 
表 1 低温 调 相 VM-PT 制冷 机 主要 结构 参数 


Table 1 Main structural parameters of VM-PT 


cryocooler with cold phase shifter 


部 件 名 称 直径 /mm 长 度 /mm 
冷 端 排出 器 24 190 
冷 端 蓄冷 器 23.6 125 

脉冲 管 蓄冷 器 外 径 18,， 内 径 8.9 85 
脉冲 管 8.5 85 
毛细 管 0.4 80 
冷气 库 外 径 51, 内 径 43 110 


2 数值 计算 模型 
e 美国 Gedeon Associates 创始 
eon 开发 的 用 于 斯 特 林 循 环 模拟 的 商用 
95 年 正式 命名 为 “Sage"， 现 在 已 经 更 新 到 
一 个 版 本 。Sage 具有 可 视 化 界面 ， 对 每 个 部 件 
独 建 模 , 各 个 部 件 之 间 通 过 质量 流 、 压 力 波 和 能 量 
流 来 进行 连接 , 可 实现 整 机 的 模拟 , 并 具有 一 定 的 优 
化 功能 。 此 外 ，Sage 计算 结果 除了 给 出 各 个 端面 的 
物理 量 外 ， 其 solution grid 文件 还 给 出 了 各 个 节点 
上 的 计算 结果 , 对 其 进行 处 理 , 可 以 得 到 各 物理 量 的 
沿 程 分 布 , 这 对 了 解 制冷 机 的 工作 原理 十 分 有 益 [7， 
本 文采 用 Sage 的 Low-T cooler Model Class( 深 低温 
模块) 对 低温 调 相 M-PT 制冷 机 进行 整 机 模拟 , 三 
个 腔 体 及 气 康 采用 Generic cylinder 模块 ,排出 右 与 
蓄冷 器 整体 采用 Constrained piston and cylinder 模 
块 , 欣 钢 器 与 连接 管道 采用 tube bundle 模块 ,室温 
狼 派 和 液 氨 档 采用 point heat source 模块 ， 蓄 冷 器 
< 用 annular canister 模块 ， 脉 冲 管 采用 compliance 


Sag 
人 Dayi 


x 人 ”tube 模块 . 计算 中 ,以 无 负荷 制冷 温度 为 目标 函数 ， 


热 腔 温度 为 300 K, 中 间 腔 温度 为 77 K. 
3 计算 结果 与 分 析 
3.1 运行 参数 优化 

数值 计算 中 ， 首 先 对 制冷 机 的 运行 参数 进行 了 
优化 计算 . 图 2 给 出 了 不 同 平均 压力 和 运行 频率 下 
制冷 机 无 负荷 制冷 温度 的 变化 规律 . 增加 平均 压力 
和 运行 频率 ， 质量 流 率 会 增加 ， 葬 冷 右 做 功能 力 增 
强 , 但 同时 蓄冷 器 损失 也 会 增加 , 因此 , 平均 压力 和 
运行 频率 均 存 在 最 佳 值 . 数值 计算 中 ， 在 运行 频率 
0.8 Hz、 平均 压力 0.9 MPa 的 运行 条 件 下 得 到 了 3.69 
K 的 无 负荷 制冷 温度 ,其 在 4.2 K 可 提供 16.24 mW 
的 制冷 量 . 值得 注意 的 是 ， 计 算 模型 中 未 计算 辐射 
漏 热 损失 , 而 Sage 软件 对 姜 冷 右 损 失 的 计算 结果 比 
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实验 数值 要 小 ,因此 实验 中 的 最 佳 运行 频率 和 最 佳 
平均 压力 一 般 要 比 计算 值 略 高 
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图 2 运行 参数 对 无 负荷 制冷 温度 的 影响 
Fig. 2 Influence of operating Parameters on no-load 
temperature 


3.2 脉冲 管制 冷 机 内 物理 量 分 布 水 
为 了 进一步 分 析 低温 调 相 VM-PT 制冷 犹 重 能 
量 流动 特性 和 耦合 机 理 ， 本 文 在 前 述 
上 , 给 出 了 脉冲 管制 冷 机 内 关键 物理 乙 疯 分 布 情况 . 
图 3 给 出 了 脉 串 管 琵 冷 器 内 气体 锋 度 的 轴 向 分 布 规 
律 。 从 图 中 看 出 ,蓄冷 器 入 国生 
度 , 但 其 冷 端 部 分 温度 梯度 镑 不 . 在 车 冷 器 中 部 , 工 
质 的 温度 已 经 低 于 6 K, 因此 在 实验 中 , 填充 在 此 温 
度 下 比 热 大 于 ErsNi 的 HoCus 有 助 于 改善 制冷 机 
性 能 . 


0.0 0.2 0.4 0.6 0.8 1.0 
菩 冷 器 相对 位 置 
图 3 脉冲 管 莹 冷 器 气体 温度 分 


Fig. 3 Distribution of gas temperature in pulse tube 


regenerator 


图 4 给 出 了 脉冲 管制 冷 机 内 时 均 声 功 、 时 均 灼 


oi 
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流 、 时 均 灶 流 的 分 布 情况 。 为 了 方便 图 形 的 绘制 ， 
将 同 轴 型 脉冲 管 示意 为 直线 结构 ， 并 忽略 了 流向 变 
化 引起 的 声 功 和 烩 六 的 符号 变化 . 大 选择 向 右 的 方 
向 为 正 ， 则 冷 端 换 热 器 出 口 至 气 库 的 声 功 和 烩 流 为 
负 号 。 

理想 状态 下 ， 声 功 只 在 划 冷 妖 和 调 相机 构 内 有 
所 耗 散 ， 在 其 他 部 件 中 没有 变化 。 从 图 4 中 的 计算 
结果 可 以 看 出 ， 声 切 依 然 主要 消耗 在 蓄冷 右 和 调 相 
机 构 内 部 , 但 其 他 部 件 内 也 由 于 流动 阻力 损失 有 少 
量 的 消耗 . 根据 热 声 理论 外， 消耗 在 蓄冷 器 内 的 声 
功 会 将 热量 不 断 的 从 冷 端 泵 送 到 热 端 ， 而 消耗 在 调 
相机 构 内 的 声 功 则 主要 起 调节 蓄冷 右 内 部 压力 波 和 
质量 流 相 位 的 作用 . 

理想 状态 下 ， 蓄 冷 历 填 料 与 工 质 之 间 无 换 热 温 
差 旦 填料 桨 给 远大 于 工 质 热 容 ,无 流 阻 损失 和 轴 向 
RCR 

答 流 并 不 为 0， 基 本 上 为 一 定 值 . 这 主要 是 

W 
的 ， 这 部 分 灼 流 将 直接 影响 制冷 机 的 制冷 量 。 由 于 
制冷 机 内 部 各 种 损失 引起 的 炉 增 , 时 均 粹 流 沿 轴 向 
不 断 增 大 并 趋 于 0. 但 在 热 端 换 热带 内 ,由 于 向 外 界 
放 热 , 炳 流 有 一 个 减 小 的 过 程 。 


一 一 时 均 声 功 
一 e 一 时 均 烩 流 


一 上 一 时 均 精 流 (WA) 


图 4 脉冲 管制 冷 机 内 能 量 流 分 布 


Fig. 4 Distribution of energy flow in pulse tube refrigerator 


图 5 给 出 了 脉冲 管制 冷 机 内 压力 波幅 值 和 体积 
流 率 幅 值 的 分 布 情况 ， 从 图 中 可 以 看 出 ， 在 进入 调 
相机 构 之 前 ， 压 力 波幅 值 基本 保持 不 变 ， 仅仅 是 由 
于 流动 阻力 引起 了 一 些 减 小 . 在 气 库 内 部 ， 压 力 波 
幅 值 约 为 0.037 MPa, 其 压 比 约 为 1.07, 远 小 于 制冷 
机 内 部 压 比 1.61, 不 会 对 制冷 机 性 能 产生 影响 . 由 于 
阻力 的 缘故 , 体积 流 率 则 沿 轴 向 不 断 减 小 . 


二 
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0 蓄冷 器 的 相位 关系 和 惯性 管 - 气 库 调 相 方式 的 相位 

| 一 一 体积 流 率 幅 值 | ”、 ”具有 很 大 的 区 别 ， 更 加 接近 于 小 孔 - 气 库 型 调 相 方 

i ] “和 式 加 .因此 ,本 文中 虽然 引入 一 段 管道 连接 脉冲 管 

] = 。 热 端 和 气 库 来 进行 调 相 ， 但 由 于 其 管 径 较 小 ， 运 行 

| 交 ”频率 较 低 ， 几 乎 没有 惯性 效应 ,只 能 起 到 阻力 元 件 
| jo% 沉 的 作用 . 加 

和 对 于 气 看 合 型 制冷 机 ， 两 级 蓄冷 器 之 间 除 了 能 

i SS 半 。 量 流 分 配 之 外 ， 还 存在 质量 流 的 分 配 ， 两 级 之 间 的 

| Noo 相位 会 相互 影响 图 6 中 冷 端 蓄冷 器 的 相位 关系 与 

兴 冰 站 天 。 热 贞 换 热 器 让 单 级 VM 制冷 机 中 的 已 有 较 大 区 别 0 这 会 加 大 制 

Em | 冷 机 调 相 的 难度 . 因此， 相位 调节 是 此 制冷 机 性 能 
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改善 的 主要 研究 方向 之 一 
图 5 “脉冲 管制 冷 机 内 压力 波 和 质量 流 分 布 改善 的 主要 研究 方 门 之 


3.4 可 用 能 损失 分 析 


Fig. 5 Distribution of pressure wave and mass fow in pulse 


tube refrigerator 参考 热 声 发 动机 驱动 的 制冷 机 Da ，VM-PT 制 
Re 
3.3 相位 分 析 4 也 


图 6 给 出 了 低温 调 相 VM-PT 制冷 机 两 级 车 冷 < 对 “ On (1 妈 ) 
器 的 压力 波 和 质量 流 之 间 的 相位 关系 图 由 于 制冷、 
机 内 部 压力 波 的 幅 值 和 相位 基本 保持 不 变 ， 因 此 和 多 
用 实 轴 来 表示 压力 波 . 图 6 中 向 量 与 实 轴 的 诡 乔 
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式 中 , 肪 表示 各 部 件 的 可 用 能 损失 ,Qn 表示 热 腔 的 
吸 热量 , Th 表示 热 腔 温 度 , Ta 为 中 间 腔 温度 。 


示 了 质量 流 和 压力 波 之 间 的 相位 关系 , , 焦 仿 穴 示 压 图 7 给 出 了 低温 调 相 VM-PT 制冷 机 中 各 主要 
力 波 的 相位 落后 质量 流 的 相位 ， 而 堪 达 内 模 表示 了 ”部 件 可 用 能 损失 的 组 成 . 图 7 中, Arice、AQw 和 AQx 


质量 流 的 幅 值 . 2 分 别 为 流动 阻力 、 不 完全 换 热 和 轴 向 导热 引起 的 可 
We 用 能 损失 . 从 图 中 可 以 看 出 , 在 各 蓄冷 器 中 ,可 用 能 
脉冲 管 蔷 冷 问 插 Dion 损失 主要 是 由 于 不 完全 换 热 引起 的 ， 其 次 为 轴 向 导 
一 并 热 引起 的 可 用 能 贺 兴 ， 而 流动 阻力 引起 的 可 用 能 损 
加 失 相对 较 小 @ 首 调 相 机 构 中， 由 于 毛细 管 管 径 较 小 ， 
气体 流 事 移 快 ， 可 用 能 损失 主要 是 由 于 流动 阻力 引 
是 05 
项 


a td xw 2 deiic 
9 ES 4o, 
1.0 区 到 4u 

实 轴 流 量 /(g/s) 
图 6 两 级 蓄冷 器 压力 波 和 质量 流 相 位 关系 


Fig. 6 Phase difference between pressure wave and mass flow 


可 用 能 损失 /名 


in two stage regenerators 


一 般 来 说 ,为 了 获得 较 高 的 蓄冷 磊 效 率 ， 压 力 


F 
[> 
> 
[> 
[> 
[> 
px 


波 和 质量 流 需 要 在 蕾 冷 器 中 部 达到 同 相 由. 但 从 图 3 

中 可 以 看 出 ， 在 脉冲 管 著 冷 器 中 部 ， 质 其 流 相 位 领 0 
先 压 力 波 接近 50 度 . 在 如 此 大 相位 角 的 情况 下 , 为 图 7 主要 部 件 可 用 能 损失 比例 

了 获得 足够 的 声 功 ， 蓄冷 器 需要 较 大 的 质量 流 » 这 Fig. 7 Exergy loss ratio of main components 


会 使 得 蓄冷 器 的 流动 阻力 损失 增 大 ， 降 低 昔 冷 器 的 
效率 。 本 文中 所 使 用 的 调 相 机 构 最 初 是 按照 惯性 管 。 3.5 制冷 量 
- 气 库 调 相 方式 设计 的 , 但 从 图 中 可 以 看 出 , 脉冲 管 图 8 根据 计算 结果 给 出 了 制冷 机 两 级 冷 头 的 制 


201705.00909v1 


chinaXiv 


694 


冷 量 曲线 . 从 图 8 中 可 以 看 出 , 一 级 冷 头 热 负荷 由 0 
增加 到 0.2 W 时 , 无 负荷 制冷 温度 由 3.69 K 上 升 到 
4.03 KK, 再 增加 至 0.5 W 时 , 无 负荷 制冷 温度 已 没有 
明显 变化 。 当 二 级 冷 头 加 一 定 热 负荷 时 ， 也 有 类 似 
的 规律 出 现 . 这 表明 ， 该 制冷 机 两 级 冷 头 在 施加 较 
小 热 负 荷 下 相互 影响 较 小 ， 可 以 在 液 所 及 液 所 温 区 
同时 提供 制冷 量 . 


17 
制冷 量 /mW 
P1 一 一 级 制冷 量 
16 上 “了 P2 一 二 级 制冷 量 
500 
网 P1 
的 14 
下 
13 
12 
3.5 4.0 4.5 5.0 5.5 
二 级 冷 头 温度 区 2 
图 8 制冷 机 制冷 量 曲 线 SS 
Fig. 8 Cooling map of the oo 
2 
4 结 论 A 


本 文 利用 Sage et 制冷 机 进行 了 数 
值 计算 研究 . 该 制冷 机 由 单 级 VM 制冷 机 气 看 合 一 
级 同 轴 型 脉冲 管制 冷 机 组 成 .为 了 减 小 来 自 室温 的 
漏 热 , 制冷 机 采用 低温 调 相 方 式 进行 调 相 . 首先 给 出 
了 各 物理 量 如 夫 流 、 炉 流 、 体 积 流 、 声 功 和 温度 的 沿 
程 分 布 , 然后 对 制冷 机 相位 关系 进行 了 分 析 , 结果 表 


明 该 制冷 机 的 相位 关系 与 惯性 管 - 气 库 调 相 方式 霸 9 


很 大 的 区 别 , 与 小 孔 - 气 库 调 相 方式 类 似 . 经 过 如 运 
行 参 数 的 优化 , 制冷 机 可 以 在 在 平均 压力 0.9 MPa、 
压 比 1.61、 运行 频率 0.8 Hz 时 可 得 3.69 K 的 无 负荷 
制冷 温度 , 在 4.2 K 时 有 16.24 mW 制冷 量 . 
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